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The complete set of two-loop renormalization group equations in general gauge field theories is 
presented. This includes the /3-functions of parameters with and without a mass dimension. 


I. INTRODUCTION 

Renormalization group equations (RGEs) provide a 
unique means in the analysis of particle physics. Com¬ 
prehensive analysis of RGEs confirmed the behavior of 
asymptotic freedom in QCD, which played a pivotal role 
in establishing a non-Abelian gauge theory for the strong 
interaction [1]. The runnings of coupling constants and 
mass parameters are crucial in global analysis of high pre¬ 
cision electroweak experiments [2]. On the other hand, 
RGEs analysis extrapolated to extremely high energy 
provides a possible, in some cases the only feasible, test 
for physics beyond the SM. The RGEs is a natural ingre¬ 
dient in the analysis of grand unification theories and 
string theories. For more than ten year, it has been 
know that gauge couplings do not unify within the SM. 
This gives extra evidence against simple grand unifica¬ 
tion theories such as SU (5) without supersymmetry, in 
addition to the non-observation of proton decay. On the 
other hand, gauge couplings seem to unify at a scale 
~ 2 x 10 16 GeV in the minimal supersymmetric standard 
model, which can be interpreted as an indirect evidence 
for supersymmetry as well as unification theories [3-5]. 

Computations of RGEs in gauge theories have been 
performed for various models to different orders of per¬ 
turbation. Persistent efforts yielded recently a four-loop 
result of the /3-function of the strong coupling constant 
[6]. Two-loop RGEs of dimensionless couplings in a gen¬ 
eral gauge theory as well as the specific case of the SM 
had been calculated long ago in a series of classic pa¬ 
pers by Machacek and Vaughn [7-9]. By introducing 
a non-propagating gauge-singlet “dummy” scalar field, 
two-loop RGEs of dimensional couplings can be readily 
inferred from dimensionless results [10, 11]. These were 
used to derive the RGEs of supersymmetric theories a 
decade later [10]. 

Recentlyfll], we re-calculated the two-loop RGEs in 
the SM, in a combination of using the general results of 
[7 9] and direct calculations from Feynman diagrams. A 
new coefficient was found in the /3-function of the quartic 
coupling and a class of gauge invariants were found to be 
absent in /3-functions of hadronic Yukawa couplings. We 
also presented the two-loop /3-function of the Higgs mass 
parameter in complete form, which provided a partial but 
useful check on the calculation of the quartic coupling. 
Whenever discrepancy with the literature appears, we 
carefully inspected relevant Feynman diagrams to ensure 
consistency. In this paper, we present the complete set 


of two-loop renormalization group equations in general 
gauge field theories. This includes the /3-functions of pa¬ 
rameters with and without a mass dimension. The results 
in [11] can then be readily reproduced. 

In section 2, we present essential notations and defini¬ 
tions, along with a discussion of the differences between 
[7-9] and our analysis. In section 3, we present the 7- 
functions of the scalar and fermion fields. In section 4, 
we present the /3-functions of dimensionless parameters 
and in section 5 those of dimensional parameters. In sec¬ 
tion 6 the results are extended to semi-simple groups. We 
conclude in section 7. 


II. NOTATIONS AND DEFINITIONS 

We start with a general renormalizable field theory 
with gauge fields V A of a compact simple group G, scalar 
fields (f> a and two-component fermion fields tpj . In section 
6, these results will be extended to semi-simple groups. 
The Lagrangian of the theory can be conveniently divided 
into three parts, 

C = Cq + C± + (gauge fixing + ghost terms), (1) 

where Cq contains no dimensional parameters and C\ in¬ 
cludes all terms with dimensional parameters. Explicitly 

7/ ( 'tpj ( 7 / 7 0 a T h.C.) — Xabcdtyatybfyc^di (2) 

where C, = ±ia- 2 - Different from [7-9], we have included 
an overall 4 factor in the Yukawa coupling terms. The 
gauge field strengths are defined to be 

Ffr = - d v V A + gf ABC V*V? (3) 

where f ABC are the structure constants of the gauge 
group and g is the gauge coupling constant. Choosing 
the standard R { gauge, the gauge field propagator is 

= 6 ab g + (1 - 0^) p. (4) 

where £ is the gauge parameter. Covariant derivatives of 
matter fields are 

D^a = d^a - ig0 A b V A (fb ( 5 ) 

Fl^j = d ^o ~ igtfkVfak ( 6 ) 
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where both 0 A b and t A k are hermitian matrices, which 
form representations of the gauge group on the scalar 
and fermion fields, respectively. Since any complex scalar 
fields can always be decomposed in terms of real ones, 
scalars in this paper are assumed to be real. 6 A are thus 
pure imaginary and antisymmetric. For late convenience, 
we define the following gauge invariants: 



b a 



c?{S) = e A c e A , S 2 (S)S AB = Tr [d A e B ] (7) 
C?{F) = t±t A bl S 2 (F)S AB = Tr [t A t B ] ( 8 ) 
C 2 {G)S AB = f A(A DjBCD (g) 


C' 2 b {R) is block diagonal for each irreducible representa¬ 
tion R( = S, F) of eigenvalue of C 2 (R). 

In this paper, we use dimensional regularization and 
the modified minimal subtraction algorithm. The renor¬ 
malized coupling constants 27 in d = 4 — 2e are related 
to corresponding bare coupling constants x k by 

°o 

x° k i-r pke = Xk + Yl < 4 n) — ( 10 ) 


where /j is an arbitrary mass scale parameter, pk = 1 ( 2 ) 
for gauge and Yukawa (scalar quartic) coupling con¬ 
stants, and a[ n ' ) are to be calculated perturbatively. The 
/3-functions of Xk are defined to be 




dxk 

dp 


e=0 


It is easy to see that 


a sr' da k ^ (i) 

Px k = 2 ^ P lXl ~d^~ ~ P ka k 


(ii) 


( 12 ) 


Perturbatively, one has 




(4tt) : 


r/4 


-B 11 


(4tt) 


(13) 


where /3 7 fc and /3 77 are one- and two-loop contributions, 
respectively. The wave function renormalization constant 
Zi of the i-th held can be expressed as 


FIG. 1: Fermion radiative corrections to scalar propagator 


The constraint on Yukawa coupling matrices imposed 
by gauge invariance is given by 

Y^a+Y^fk+tfrY,%=0. (17) 

When a fermion is involved simultaneously with a 
Yukawa coupling and a gauge coupling, one usually has 
the combination Y^tfj. + t A !*Y lk . This is to say, in these 
cases, Y a is always preceded by a t A * or followed by a 
t A . Eq. (17) can thus be used to simplify the gauge 
structure of Feynman diagrams. In [7-9], fermion fields 
are implicitly assumed to be real so t A are pure imag¬ 
inary and antisymmetric. t A * was substituted by —t^ 
and eq. (17) is reduced to 

Y b X a +Y° l t A -t A Y l l = 0. (18) 

However, in most physical theories, fermions are in gen¬ 
eral Weyl and complex, with only possible exceptions of 
neutrinos and gauginos as real Majorona fermions. For¬ 
tunately, close inspections of Feynman diagrams show 
that both approaches yield the same final results in most 
cases. In cases where they differ, we will retain t A *. 

By assuming real fermion fields, the direction of a 
fermion propagator needs not to be discriminated either. 
Shown in Fig. 1 are two fermion loop diagrams which 
contribute to the propagators of scalar fields. The two 
diagrams are different for complex fermions and the total 
result is proportional to 

Yf(S) = ^Tr (Y +a Y b + Y +b Y a ). (19) 




(14) 


The corresponding anomalous dimension is 

Id 1 ^ ■* 

7i = 0^, X ° gZi = 


2 ^ dp 

Also perturbatively, one has 


i 


dxi 


1 




.ii 


(4tt) 


2 ii 


(4tt) 


4 /1 


(15) 


(16) 


where 7 / and 7/ 7 are one- and two-loop contributions, 
respectively. 


Y 2 b (S) forms a hermitian matrices and is block diagonal 
for each irreducible representation of eigenvalue Y 2 (S). 
For real fermions, the two terms in Y 2 ab (S') become equal 
and they reduce to a common factor Tr(Y + °Y b ). Again, 
we will retain eq. (19) and other similar combinations. 
The second part of the Lagrangian is, 

1 7~n? 

2 “I” h.C.] — (fiafib 

-^4>aMc (20) 

In principle, the (3 functions of (nif)jk, m 2 ab and h abc can 
be calculated directly. But this would be tedious and can 
be avoided. To do so, we introduce a non-propagating 
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dummy real scalar field q with no gauge interactions 
and rewrite the Lagrangian as, 

A = +/i-c) - 

( 2 !) 

with the substitution of Yf- = (rrif)ij, A ab jj = 2 rn 2 b and 
A abed — habc- The (3 functions of m/, ?n 2 b and h abc are 
equal to those of the new Yukawa coupling Y d , quartic 
scalar coupling A js, and A abcd ~. The latter can be read¬ 
ily obtained from the (3 functions of Y a and A a bcd by 
suppressing both the summations of the d-type indices 
and related gauge couplings. 


III. WAVE FUNCTION RENORMALIZATION 
A. Scalar wave function renormalization 


To one loop, the anomalous dimensions of scalars are 
lab = 2 nY 2 ab (S) - g \3 - £)Cf(S) ( 22 ) 

To two loops, 


lib = -g Cf(S) 


f-2 Z-\e)Ci(G) 


™kS 2 {F) ^ ^S 2 (S) 




+ - g*C$°(S)C?(S) - 3 K H 2 ab (S) 
- 2 K H 2 ab (S) + lO K g 2 Y* b F (S) 


(23) 


where k = 1 / 2 ( 1 ) for two(four)-component fermions, and 

A„j(iS) = -Xacde^bcde (24) 

0 

H Ib(S) = ^TiiY a Y +b Y c Y +c + Y +a Y b Y +c Y c )(25) 
Hl b {S) = ^TiiY a Y +c Y b Y +c + Y +a Y c Y +b Y c )(26) 
Y 2 b F (S) = ^Tr [C 2 (F)(Y a Y +b + Y b Y +a )\ (27) 

In contrast with [7], Y“ b (S), H% b (S), H% b (S) and Y 0 2 b F (S') 
are modified so they are hermitian for complex fermions. 
As mentioned above, it is due to that complex fermion 
lines have two distinct directions in Feynman diagrams. 
The contributions of the two parts are hermitian conju¬ 
gate to each other. For real fermions, both contributions 
are real. So they are equal and can be combined. Similar 
structure appears in later results. Here and hereafter, we 
express the Casmir factors as C 2 b (S) instead of C 2 (S)S ab 
and C 2 oc (S)C'| b (S) instead of C 2 (S)S ab , to accommodate 
reducible representations. 


B. Fermion wave function renormalization 


To one loop, the anomalous dimensions of fermions are 


if = \y a Y +a + g 2 C 2 (F)^ 

To two loops, 

7 fj = -}-Y a Y +b Y b Y +a - lnY a Y +b Y£ b (S) 
8 2 

2 ^ Cf{S)Y a Y +b - jC 2 (F)Y a Y +a 


(28) 


9 


1 , 


- -Y a C 2 (F)Y +a 


9 A C 2 {F) 


f+ 2 t+\e)C 2 {G) 


1 


- 2 kS 2 (F) - -S 2 (S) 


-| 9 4 [C 2 (F)] 2 (29) 


Here the only modification is in the Y 2 ab (S'), which ap¬ 
pears in the first line of eq.(29). This is due to one-loop 
fermion sub-diagrams in scalar field propagators. 


IV. DIMENSIONLESS PARAMETERS 


A. The gauge coupling constant 


The /3-function for the gauge coupling constants has 
been extensively studied. Recent results are up to the 
fourth order. Here we include the two loop results for 
completeness, 


P(g) = - 




( 47 t ) 2 
2k 


yG 2 (G) - | kS 2 (F ) - ±S 2 (S) 


+ (M^> > - 


r 


( 47 t ) 4 


34 o 

y[C 2 (G)] 2 


— K 


20 

4G 2 (F) + —G 2 (G) 


2G 2 (5) + -G 2 (G) 


where Y±(F) is defined through 

1 


S 2 (F) 
S 2 (S) 


M p ) = yy Tr [C 2 (F)Y a Y+ a ] 


(30) 


(31) 


d{G) is the dimension of the group. Note that though 
the Yukawa couplings are normalized differently, eq. (30) 
assumes the same form as eq (6.1) of [7]. 


B. The Yukawa couplings 

The /3-functions of the Yukawa couplings can be ex¬ 
pressed as 

/3 a = 7 a + 7 +F Y a + Y a 7 F + 7 f b Y b (32) 
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where 7 0 are the anomalous dimensions of the operators 
< / > aVh'C' i /’fc> 1 F an d lab are anomalous dimensions of 
the corresponding fermions and bosons, respectively. To 
one loop, 


0? = i [Y+(F)Y a + Y a Y 2 (F)\ +2Y b Y +a Y b 

rb\rab 


+2 K Y b Y 2 ab (S) - 3g 2 {C 2 (F), Y a } 


where 


Y 2 {F) = Y +a Y a 


(33) 


(34) 


To two loops, 


fill = 2Y c Y^Y a (Y' JrC Y* ) — 

- Y h [F 2 (F)y +a + Y +a Y+{F)] Y b 

- - [Y b Y 2 {F)Y +b Y a + Y a Y +b Y+{F)Y b ] 

8 

- 4 nY 2 ac {S)Y b Y +c Y b - 2 nY b H 2 b (S) 

- ^nY 2 bc (S)(Y b Y +c Y a +Y a Y +c Y b ) 

- 3 K Y b H 2 ab (S) - 2\ abcd Y b Y +c Y d + l -A 2 ab {S)Y b 

+ 3 g 2 {C 2 (F),Y b Y +a Y b } + 5 g 2 Y b {C 2 {F), Y +a }Y b 

- 7 -g 2 [C 2 (F)Y+(F)Y a + Y a Y 2 (F)C 2 (F)} 

- jg 2 [Y b C 2 {F)Y +b Y a + Y a Y +b C 2 {F)Y b } 

+ 6 g 2 H a 2t + 10 Kg 2 Y h Y 2F {S) 

+ 6 g 2 [C b 2 c {S)Y b Y +a Y c - 2 C% c (S)Y b Y +c Y b ] 

+ ^g 2 C% c (S)(Y b Y +c Y a + y a Y +c Y b ) 

- \g\[C2{F)] 2 ,Y a } + 6 g 4 Cf(S){C 2 (F),Y b } 

+ 5 4 -yC 2 (G) + y«S 2 (F) 


11 


S 2 (S) 


12 

g 4 C% b (S) 


{C 2 (F),Y a } - ^-g 4 C? (S)C bc (S)Y c 


40 1 

-C 2 (G)-2kS 2 (F)--S 2 (S) 


where 


= t A *Y a Y +b t A *Y b + y b t A Y +b Y a t A 


Y(3 5) 


(36) 


The second term in the fourth line, the first term in 
the sixth line, and the second term in the tenth line, 
are expressed in terms of the liermitian gauge invariants 
H 2 b (S), H 2 b (S), and Y 2 b F (S), respectively. They are all 
from the 7 -functions of the scalar fields. They differ from 
the fourth line and the ninth line of (3.3) in [ 8 ], which 
included only parts of the expressions. In H 2t , in front 
of Y a and Y b one has t A * instead of t A , again due to the 
fact that fermions are complex. 


C. The scalar quartic couplings 

The /3-functions of the scalar quartic couplings can be 
expressed as 


Pabcd — 7 abed T ^ v 7 ( F) Aabcd 


(37) 


where labed are the anomalous dimensions of the opera¬ 
tors (f)a 4 *b$c 4 *d,i 1 S {i) is the anomalous dimension of the 
scalar field i. To one-loop, 


Pabcd Aabed I abc< j T 2h'A abcd 

—3g A abcd + 3 g A abcd 


where 


A ,abed — o ^ 1 A abe fA, 


efed 


H abcd = Tr (y a Y +b Y c Y +d ) 

perms 

A a bcd = Y 2 (j) Agbcd 

i 

— y ^ G 2 [PjAgbcd 


\S 

^abed 


A abcd = ^ E {0 A ,0 B }ab{e A ,e B } cd 


To two loops, 

Pabcd = y ^ P) Aabcd ~ A abed ~ 4ft A a 


(38) 

(39) 

(40) 

(41) 

(42) 

(43) 


2 Y 
bed 


« \ § H* bcd - E t 3i?2 (*) + 2 H 2 (i)\ A a bcd 


+ 4 K(H^ bcd + 2 H^bcd + 23?o bed) 

+ g 2 {2A 2 a s bcd -6A 2 a l cd + MH B 


_ ttF \ 

abed ^abed) 


+10KY / Y 2F (i)X abcd 


g 


Y c 2(G)-y kS ^f)-^s 2 (S) 

-~A x 

Q^^abcd 2 “^abed abed 

+^ K (Bl bcd — 10 Flbcd)} 

6 j ^G 2(G ) - yftS 2 (F) 7 -S 2 (S ) 


1 

A s _i_ 07 a9 

2 ^abed ' ^ ‘ ^abed 


\ s 
^abed 


-A-abcd 

(44) 


where A 2 (i) 1 H 2 (i ), H 2 (i) and Y 2F (i) are the eigenvalues 
of the invariants of A 2 ab {S), H 2 b (S), H 2 b (S), Y 2F (S), 
respectively. The other invariants are defined by 


A a bcd — 4 y ] A abe fA ceg hA d fgh 


(45) 


perms 
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\ 2Y 
abed 

= l E Y 2 f9 (S)\ abef \ cdeg 

(46) 


perms 


ttX 

abed 

= l E KbefTr(Y c Y +e Y d Y+f 



perms 



+ Y Ac Y e Y +d Y?) 

(47) 

ttY 

f^abed 

= E Tr [Y 2 (F)Y +a Y b Y +c Y d ] 

(48) 


perms 


ttY 

abed 

_ _ y ^ 



^ perms 



_|_y+eyay+ey6y+Cyd j 

(49) 

rr3 

abed 

_ _ y ^ y r (y a y+feyey+cydy+e) 

(50) 


perms 


a 2S 

abed 

= l E Cl 9 (S)XabefXc d e g 

(51) 


perms 


A 29 

^abed 

= g 'y ' y X abe f X cdgb 0 eg 9f h 

(52) 


perms 


ttS 
**abed 

= E F! 2 {i)H abcd 

(53) 

ttF 

abed 

= E Tr [{C 2 (F),Y a }Y +b Y c Y +d ] 

(54) 


perms 


A ss 
abed 

= E 1^2 (^’)] 2 X abcd 

i 

(55) 

A x 

-™-abcd 

= \ E A «be f {e A ,e B } ef {o A ,e B } cd 

(56) 


perms 


A x 

abed 

= \ E *a b ef{0 A ,0 B } Ce {0 A ,O B } df 

(57) 


perms 


tdY 

F* abed 

= i E( 0 A ,0 B }abYr[t A *t B *Y c Y +d 



perms 



+Y c t A t B Y +d ] 

(58) 

by 

F*abed 

= J E {0 A ,0 B }abTr{t A *Y c t B Y +d ) 

(59) 


perms 


A s 

^abed 

= ^ ^ C2 (fyAgbcd 

(60) 


i 


K bcd = \f ACE f BDE E {^ B U{&°,B D U 61) 

perms 


In the first term of B^ bcd , Y c is preceded by t A *t B * in¬ 
stead of t A t B , since the fermions are complex. Similarly 
in B^ bcd , Y c is preceded by t A * instead of t A . Since there 
is only one t factor in this case, this introduces one extra 
minus sign. Therefore, the relative sign between B^ bcd 
and B^ bcd in Eq.(44) is minus while it was plus in [9]. In 
addition to H 2 b (S), H 2 ab {S ), and Y 2 b F (S), H^ bcd is also 
re-expressed to be hermitian. 


V. DIMENSIONAL PARAMETERS 
A. Fermion mass 

The /3-functions of fermion mass can be inferred from 
those of the Yukawa couplings by taking the a-indices to 
be dummy. The trilinear scalar terms start to contribute 
from two-loops. The one loop result is 

PL f = \ [Y 2 + (F)m f + m f Y 2 (F)\ + 2Y b m+Y b 

+KY b Tr(m+Y b + m f Y +b ) 

-3 g 2 {C 2 (F),m f } (62) 

The two loop result is 

(3 I J Lf = 2 Y c Y +b m f (Y +c Y b -Y +b Y c ) 

- Y b W 2 {F)m+ + m+Y+(F)] Y b 

- - [Y b Y 2 (F)Y +b m f +m f Y +b Y+(F)Y b ] 

8 

- 2KY b Y +c Y b Tr(mjY c + m f Y +c ) 

- ^KY£ c {S)(Y b Y +c m f +m f Y +c Y b ) 

- -nY b TV \y 2 (F)Y +b m f + mtY+{F)Y b 

2 L J 

- KY b Tr(Y c m^Y c Y +b + Y +c m f Y +c Y b ) 

- 2 h bcd Y b Y +c Y d + -h cde X bcde Y b 

+ 3 g 2 {C 2 (F),Y b m+Y b } + 5 g 2 Y b {C 2 (F), m+}Y b 

- - A g 2 [C 2 (F)Y+{F)m f + m f Y 2 (F)C 2 (F)\ 

- ^g 2 [Y b C 2 {F)Y +b m f + m f Y +b C 2 (F)Y b } 

+ Qg 2 [t A *m f Y +b t A *Y b + Y b t A Y +b m f t A } 

+ 5ng 2 Y b Tr[C 2 (F) (m f Y +b + Y b mf)] 

+ 6 g 2 C bc (S)Y b m+Y c ~ \g A {[C 2 {F)f , m f } 

+ | g 2 C b 2 c {S){Y b Y +c m f +m f Y +c Y b ) 

+ g 4 ——C 2 {G) + —kS 2 (F) 

+ §S 2 (S) {C 2 (F),m f } (63) 

B. Trilinear scalar couplings 

The /3-functions of trilinear scalar couplings can be in¬ 
ferred from those of the quartic couplings by taking one 
of the four indices to be dummy. The fermion masses 
contribute from one-loop. One loop result is 

PLc = A lbc - 8KH abc + 2nK y abc - 3g 2 A B bc (64) 
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where the invariants are defined as 


A 


— _ ^ A abefh, 


abc 


efc 


(65) 


perms 


H abc = i E Tr (m f Y +a Y b Y +c + Y a m+Y b Y+{$6) 


2 Y 
abc 


A,abc — (i)habc 

i 

A, a b c = ^ ' C 2 (l)h abc 

i 

The two loop result is 

#1 = ^E a 2 (^-a 3 o 6 c - 4 k a: 

i 

+ Kjs^ + s^V 

~E[3ff 2 «+2ff 2 (i)] /i a6 c| 

+ ^ K (H Y bc + 2 ff^, c + 2H^ bc ) 

+ .9“ |2A^f c - 6 A 2 ^ c + 4 K(H B bc — H B bc ) 

+io K^y 2F (i)ii*| 


(67) 

( 68 ) 


5 


4 |J|c 2 (G)-y^ 2 (F)-^ 2 (5) 

_£a ss --4 a --A x 

2 lv abc 2 ^abc abc 

+4«(BL - lO^L)} 

where the invariants are defined as 
1 


A: 


be 


(69) 


^abc — 9 E/ (^■abefA. ce glhfgl+X aeg lXbfglh ce f\7 0) 


perms 


A 2 L = ^ E Yi 9 {S)Xa bef h c 


(71) 


perms 


Habc = o E A obe /Tr (F c F+ e m / y+/ 


perms 


+m f Y +e Y c Y +f + y +c F e m^y / 

+m^y e y +c y / ) 

H h abc = J E /i ae/ Tr (y b y +e y c y +/ 

perms 

_j_y+6yey+cy / \ 

tf/L = E Tr E(F)(m+y a y +b y c 


(72) 


(73) 


+y +0 TO/y +b y c + Y +a Y b mp c 
+Y +a Y b Y +c m f )\ (74) 

tf/L = \ E Tr (y e m+y e y +a y & y +c 

perms 



+Y e Y +a Y e mp b Y +c 

+Y e Y +a Y e Y +b m f Y +c 

+ yey+ayey+byc TO + + ^ 

(75) 

tt3 

n abc 

= i E Tr(?u / y +a y e y +b y c y +e 



perms 



+Y a m.p re Y +b Y c Y +e 
+Y a Y +b Y e mp c Y +e 
+ Y a Y +b Y e Y +c m f Y +e ) 

(76) 

x 2S 

abc 

= \ E Cl 9 (S)h aef X bceg 

(77) 


perms 


A 2 # 

^abc 

= \ E h aefhcgl0? g 0fl 

(78) 


perms 


ttS 

^abc 

= E C 2 (*)#a6c 

(79) 

ttF 

u abc 

= E Tr[{G 2 (F),m/}y +a y b y +c 



perms 



+{C 2 (F),Y a }mp b Y +c 
+{C 2 (F),Y a }Y+ b m f Y +c 
+ {C 2 (F),Y a }Y+ b Y c m f ] 

(80) 

\SS 

abc 

= Kbc 

i 

(81) 

■™-abc 

= \ E Ke f {o A ,e B } ef {9 A ,e B }b C 

(82) 


perms 


■™-abc 

= \ E h »ef{0 A ,0 B } be{0 A ,0 B } cf 

(83) 


perms 


n Y 

abc 

= \ E {0 A ,0 B } a bTr(t A *t B *m f Y +c 



perms 



+m f t A t B Y +c + t A *t B *Y c m+ 

+Y c t A t B m 

(84) 

f>Y 

abc 

= l E {0 A ,0 B }abTr(t A *mft B Y +c 



perms 


+t A *Y c t B mf s ) (85) 


C. Scalar mass 


The /3-functions of scalar masses can also be inferred 
from those of the quartic couplings by taking two of the 
four indices to be dummy. From one loop, both fermion 
masses and trilinear terms contribute. One loop result is 

drn\ ~ ^ef Aabef + haefhbef AhvH ab 

-3g 2 A s ab + 2nA Y ab (86) 
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where the invariants are defined as 

H ab = Tr[ (Y a Y +b + Y b Y +a )m f mj 
+(Y +a Y b + Y +b Y a )m+m. f 
+Y a m^Y b m + f + m f Y +a m f Y +h (87) 

A ab = ( 8S ) 

i 

a Y ab = (89) 

i 

The two-loop contribution is 
0m* ab = ^Y^ A2 ^ m ab- \ A lb~^ A2 ab 


+ k ^4H x b - £ i 3H2 W + ^W] m2 b 

+ ^ K (H^ b + 2 H% b + 2 H 3 b ) 

+ g 2 {2Af b - 6A % + 2 n{H s ab - H B ) 


+ 10kY y2F ^ m lbj 

- ff 4 { f G 2 (G) - y«5 2 (F) - ^S 2 (S) A s ab 


--\ ss -^A x - -A x 
2 ah 2^ ab 

+2 K (Bl b - 10 B^)} 


where the invariants are defined as 


A ab — X a bef heglh f g l -f“ 2lln e f\ ae gl\bfgl 
+2 h 

ae fh/giXbegi + 2,h be fhfgi\ aeg i (91) 

A ab ~ ^ 2 9 (S)(m eg \ a bef + haefhbeg) (92) 

= \^.W’n,Y+’m,Y+’ + h.c.) 

+m 2 ef TiiY a Y +e Y b Y +f + h.c.) 
+h ae fYr(Y b Y +e m.fY + f + h.c.) 
+h bef Tv(Y a Y +e m f Y +f + h.c.) (93) 

HU = 2Tr \ i {Y 2 (F),mjm,f}(Y +a Y b + Y +b Y a ) 

+ 2Tr \Y 2 (F)Y +a mf(Y +b mf + m.jY b ) 
+Y 2 (F)m^Y a {Y +b m. f + m^Y b ) 
+Y 2 (F)Y +b m f {Y +a m f + rrijY a ) 
+Y 2 (F)mjY b (Y +a m f + mjY a ) j (94) 

Hll = Tr (Y e Y +a Y e Y +b + Y e Y +b Y e Y +a )m f m^ 
+Y e m^Y e m^{Y a Y +b + Y b Y +a ) 
+Y e Y +a Y e mj(Y b m + + m f Y +b ) 
+Y e m~fY e Y +a (Y b mj + m f Y +b ) 


+Y e Y +b Y e m+(Y a m+ + m f Y +a ) 
+Y e mfY e Y +b (Y a mf + m f Y +a ) + h.c. 1(95) 

Hab = Tr \{Y a Y +b + Y b Y +a )Y e m^m f Y +e 
+m f m^Y e (Y +a Y b + Y +b Y a )Y +e 
+Y a m+Y e (Y +b m f + m+Y b )Y +e 
+m f Y +a Y e (Y +b m f + m~j Y b )Y+ e 
+Y b m^Y e (Y +a m f + mjY a )Y +e 
+m f Y +b Y e {Y +a m f + m^Y a )Y +e 1 (96) 

A ll = Cl 9 (S)\ abef m 2 eg + Cl 9 (S)h aef h beg (97) 

A l 9 b = X abef m 2 gl + haefhbgiOtgOf,, (98) 

Hab = Y° 2 ^ H -b (99) 

i 

Hab = 2Tr[ {C 2 (F),Y a }Y+ b m f m+ 
+{C 2 {F),Y b }Y+ a m f m+ 
+{C 2 (F),m f }m+(Y a Y+ b + Y b Y +a ) 
+{C 2 (F), Y a }m+(Y b m+ + m f Y +b ) 
+{C 2 (F),m f }Y +a {Y b m+ +m f Y +b ) 
+{C 2 (F), F b }rn+(F a TO+ + m f Y +a ) 
+{C 2 (F), m f }Y +b {Y a m + f + m f Y +a ) 1 (100) 

Af b = ^ [CM*)] 2 ml,, (101) 

i 

A ab = m 2 ef {0 A ,e B } ef {9 A ,8 B }a b (102) 

A x ab = m 2 ef {e A ,e B } ae {0 A ,0 B } bf (103) 

B Y ab = {e A ,e B } ab Tr(t A *t B *m f mj + m f t A t B mp.04:) 

Bib = {0 A ,0 B } ab Tr(t A *mft B mf) (105) 

VI. EXTENSION TO NON-SIMPLE GROUPS 

So far, the gague group is assumed to be simple. These 
results can be extended to semi-simple groups by assign¬ 
ing appropriate substitution rules, based upon close in¬ 
spections of relevant Feynman diagrams [7-9]. 

Assume the gauge group is a direct product simple 
groups, Gi x • • • x G n , with corresponding gauge coupling 
constants gi,... ,g n . The substitution rules for the gauge 
coupling constants are 


g 3 C 2 (G) g 3 k C 2 (G k ) (106) 

g 3 S 2 (R) -> g 3 S k (R) (107) 

g 5 [C 2 (G)} 2 -> gl[C 2 {G k )) 2 (108) 

g 5 C 2 (G)S 2 (R) —> g 5 k C 2 (G k )S k (R) (109) 


g 5 C 2 (R)S 2 (R) -> Yti9?C l 2 (R)S k (R) (110) 

l 



here and hereafter, k and l are subgroup indices, R can 
be either S or F. 

For other (3 and 7 functions, we first have the following 
general substitution rules 


g 1 2 3 C 2 (R) 

- 

Yal^iR) 

k 

(in) 

g 4 5 C 2 (G)C 2 (R) 

- 

Ygtc 2 (G k )c k (R) 

k 

( 112 ) 

g 4 S 2 (R)C 2 (R') 

- 

Y9tS k 2 (R)C k (R') 

k 

(113) 

g 4 C2(R)C2(R') 

- 

Y9 2 k9?C k (R)C l 2 (R') 

k,l 

(114) 

In HZt, Kh cch in¬ 
stitution rules are 

B ab > B abcd’ B ahc> B ab > * h e Sub- 

e A -> 

9t 

e B -> ef 


t A - 

At 

t B - tf 


(it A * -> 

+A* 
l k > 

t B * -> tf*) 
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4 

2 2 

9 k 9i 

(115) 

For example, g 4 B\ bcd is substituted by 



g 6 S 2 (R)A abcd by 

-J2ti9?S k 2 {R) Y. (A k ac,efKfM+ A ae,fdKb,cf) 

k,l perms 

(121) 

g 6 C 2 (G)A abcd by 

4 55 9k9lG 2 (Gk) 55 (Aac,e/Ai/,6d + A k e j d A l ebc f) 
k,l perms 

( 122 ) 

and g 6 A s abcd by 

Y 55 SfeCf (*) J Y 9 ? 9 m 

k i l,m 

X Y ( A ac,efKfM + <e,fdK b ,cf) (123) 
perms 

Finally, one has 

9 6A9 abcd -Y3 6 AcAk) ( 124 ) 


\Yti9l Y {OktfUTr(t A *tf*Y c Y+ d +Y c t A tfY+ d ) 

k,l perms 

( 116 ) 

In quartic coupling, trilinear coupling and scalar mass- 
square terms, further substitution rules are needed. In¬ 
troduce a new tensor 

A ab,cd = 0 0 A ) ac (0 A ) bd (117) 

so the gauge invariants A a f, C( i can be rewritten as 

Aabcd — ^ ^ ^ (-^-ac,e/ ■Aef,bd A A ae ^f d A ebc f). (H 3 ) 

perms 

The substitution rule for A ab , c d is 


VII. CONCLUSIONS 

We have now presented the complete set of two-loop 
renormalization group equations in general gauge theo¬ 
ries. This included the /3-functions of parameters with 
and without a mass dimension. We have so far restricted 
the gauge groups to be senri-sinrple. If the gauge group 
contains more than one U( 1) group, the situation is sub¬ 
tle. If there are no kinetic mixings between the U( 1) 
groups after renormalization, these results can applied 
by straightforward extension. If there are kinetic mix¬ 
ings, these results have to be applied with caution. In 
general, modifications will be warranted [12], In the case 
of the SM where one has only one U{ 1) group, one readily 
reproduces the results in [ 11 ]. 


g A ab ,cd 55 9k A ab,cd (H9) 

k 
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